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7.1 Introduction 

L’empilement de plusieurs plis de differentes orientations qui sont colles ensemble forme 

un stratifie. 


Definition du code de description 

Chaque pli est note par 1’ angle (en degre) entre la direction longitudinale des fibres 
et l’axe x. 

Des plis adjacents sont separes par une barre oblique (/) si leurs angles ont de 
differentes valeurs absolues. 

La sequence d’empilement commence a partir du premier pli superieur et se termine 
au dernier pli. 

Un indice est utilise afin d’indiquer le nombre de plis adjacents possedant la meme 
orientation. 
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Figure 7.1 Quelques exemples d’empilement 
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Definitions 

Stratifie symetrique 


Stratifie equilibre 


Stratifie orthogonal 


: Les plis sont disposes symetriquement par 

rapport au plan moyen geometrique (plan (x,y) 
passant par le milieu de l’epaisseur du stratifie). 

: Nombre de plis a +9 est egal au nombre de 

plis a -0 


Plis de 0° et de 90° 


z 
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1.5 


Stratifie non equilibre 
et non symetrique 



z Stratifie non equilibre 
et symetrique 


-1.5 

-0.5 
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Stratifie equilibre 
non symetrique 



-0.5 

0 
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symetrique 


-1.5 
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Figure 7.2 Differents types de stratifie 


Notes : 

• Un pli dont les axes naturels ne coincident pas avec les axes de reference se 
comporte comme un stratifie non equilibre mais symetrique. 

• Pour qu’un stratifie soit symetrique et equilibre, il doit comprendre au moins 
quatre plis. 

• Un stratifie symetrique et equilibre doit comprendre un nombre pair de plis. 
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7.2 Theorie (classique) des plaques stratifiees (avec couplage) 

Plusieurs plis a orientation quelconque. 

Empilage des plis symetrique ou non. 



Figure 7.3 Systeme de coordonnees et les elements de reduction [Gibson] 


(x, y) est le plan moyen geometrique (PMG) 

u,v etw sont les deplacements suivant les directions x, y et z. 

Hypotheses de base : (Etat plan de contraintes et petites deformations ) 

Element est constitue par un nombre arbitraire de plis parfaitement lies les uns 
aux autres et dont les axes naturels ne coincident pas obligatoirement avec les axes de 
references. 

L’epaisseur t de l’element est tres faible par rapport aux deux autres dimensions 
(longueur et largeur). 

Les deplacements u, v, w sont faibles par rapport a l’epaisseur t de l’element. 

Les deformations e x , £ y et y xy sont faibles. 

Les deformations de cisaillement et y yz sont negligeables. 

Les deplacements u et v sont des fonctions lineaires de z. 

La deformation normale £ z est negligeable. 

Les relations a - s dans chacun des plis sont conformes a la loi de Hooke. 
L’epaisseur t de l’element est constante. 

Xxz et x yz sont nulles a la surface de l’element (z = ± — ). 
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Elements de reduction 

Contraintes resultantes 


Considerons un stratifie constitue de 3 plis d'une epaisseur de 1 mm chacun : 


-<—£ = 10 mm— ► 

1 -1.5 


1 mm 

acier 

-0.5 

1 mm 

aluminium 

+0.5 

+1.5 

1 mm 

acier 

t 






z 



e y = 0.07 % 


4 

mu 


A 


► 

A 

'///////////A 




Figure 7.4 Elements de reduction 


Equilibre : 2F ac + F al = P y 


A^ac * -^ac) (Aal * -^-al) A (3) 


Compatibilite : s ac = s al (2) 


Rel. a - s : 


< 


E, 



( 3 ) 


(3) dans (2) => ^ 
E_ 


= a. 



70GPa - p ac (4) 
210 GPa 3 
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(4) dans (1) => 2a ac x A ac +^x A al 


or A ac = A al =10 mm x 1mm = 10 mm 2 


d'oii 



on obtient : 


a ac =150 MPa 
a al = 50 MPa 


et 


s =^xl00= 150 xlOO = 0.07% 

y E al 210xl0 3 


g x n'est pas constante a travers l'epaisseur comme dans les materiaux isotropes. II faut 
choisir une autre grandeur pour caracteriser l'etat de charge du stratifie. 

3 

Equilibre des forces : P y =Z( a y)k A k 

k=l 

or A k = t k x l k etl k =c te =l 


d'ou : 


P y =Z( a y)k xt k xl k = 1 Z( a y)k xt k 


k=l 


k=l 


On definie la contrainte resultante suivant la direction y comme : 


N y =-f = Z( a y)k xt k 

I 1,_1 


P„ 2 


k=l 


N y est done 1’ element de reduction de contrainte o y dans l'epaisseur 


n 


N y = Z(°>)k xt k 


k=l 


(a y )k peut etre une fonction de z comme dans le cas de flexion. 
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N y represente la charge totale suivant la direction y, supportee par un stratifie d'epaisseur 

f r- ivt A 


t et de largeur unitaire dans la direction x. 


Force 


N 


unrte = 
v Largeur mm 


= MPa - mm 


y 


De faqon plus generale, N y s'ecrit : 


N y =Pt CT y dz = L{C_ 1 (c 'y ) k dz | 

2 k — 1 

(7.1) 

1 N f y 

N x = Pt a x dz = £ {P 1 (°x)k dz j 

2 k— 1 

( 7 - 2 ) 

1 N / \ 

N xy = Pt °xy dz = £ (PP ( a xy)k ( dz }} 

(7.3) 


Les moments resultants de flexion et de torsion s'ecrivent de la meme faqon : 

M x =Pt a x zdz= Zlf Zk ( a x)k xzdz j ( 7 -4) 

~2 k=l l Z k-1 J 

M y = P t a y zdz = Z jf Zk ( a y)k xzdz l ( 7 - 5 ) 

-- k=l l Z k-1 J 

M xy =Pt T xy zdz= zl lp ( X xy)k xzdz l ( 7 -6) 

~2 k=l l k-i J 

T , . , , , , n . N - mm 

L unrte de mesure My est (Force x hauteur / largeur) = = N 

mm 

Les trois Ny et trois My forment un systeme statiquement equivalent de l'etat de 
contrainte du stratifie, mais qui est applique au plan moyen geometrique. 


Relations entre deformations et deplacements 

Les deplacements du plan moyen geometrique seront designes par u° et v° suivant les 
directions x et y, et par w suivant la direction z. 
u° = u°(x, y) 

v° = v°(x, y) (7.7) 

w = w(x, y) 
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L’hypothese des faibles deformations permet d’ exprimer les deformations du plan moyen 
geometrique : 


1 


du 0 

dx 

dv 0 

dy 

du n 


(7.8) 


xy 


+ 


dv n 


dy dx 


Le probleme consiste a determiner les deformations en un point p situe a une distance 
quelconque z du plan moyen geometrique 



Figure 7.5 Section droite d’un stratifie avant et apres Faction de la charge 
exterieure [Hyer] 

Sous Faction de la charge externe, le point p° qui est situe dans le PMG se deplace d’une 
distance u° suivant la direction x. Selon l’hypothese de Love-Kirchoff (les sections planes 
restent planes), la ligne pp° reste droite et perpendiculaire au PMG. Le point p° se deplace 
d’une distance : 


Or 


U=U 0- Z /?x 

tgPx ~ (3* 


<9w 


dx 


Sw 

done : u = u n - z 

0 dx 

De la meme faqon, on obtient le deplacement du point p suivant la direction y : 

Sw 


v = v 0 -z 


dy 


(7.9) 

(7.10) 
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En utilisant les relations entre deformations et deplacements, les composantes de la 
deformation en un point p quelconque de l’epaisseur ont pour expression : 


du d , o <3w\ du L 


z x = — = —(u -z— ) = 


8x ox 


OX dx 


-z- 


d 2 w 


dx 


2 = £ x~ Z 


d 2 w 


dx' 


(7.11) 


dv 5, 0 5w\ dv u 




dy dy 


— z 


d 2 w 


d 2 w 


dy dy dy 


2 =s y _z 


dy 


0 T 

Yxv=Yxv- 2z 


I xy i xy 


d 2 w 
dxdy 


(7.12) 

(7.13) 


Le rayon de courbure du plan moyen geometrique peut s’ exprimer en fonction de w : 
1 d 2 w 


= k„ 


p dx^ 

d’oii :s„ =s° + zk v 


(7.14) 


de la meme faqon on obtient : 


£ y - £ y +ZK y 


on: K y =- 


d 2 w 
dy 2 


et v = Y + zk 

< xy / xy xy 

„ d 2 w 

ou : k = -2 

xy dxdy 

Le champ des deformations en un point quelconque de l’epaisseur d’un stratifie peut 
s’ecrire sous forme matricielle : 



r > 


'8°' 


f 1 






K x 


< 

8 

r = i 

s u 

f + z 1 

K 

> 


y 


y 


y 



[YxyJ 


r 

§ ° 
v 


KyJ 



(7.15) 
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Relation entre les elements de reduction et les deformations 


( 1 ) 

\ 2 l 


/ 3 ^ 



/ Plan moyen geometrique \ 

^3 

Z2 

*1 


r 

t_ / 


Z / 



Zk 

Zn- 1 

2 

k 










N 



Figure 7.6 Geometrie et cotation des plis d’un stratifie [Gibson] 

Dans le k e pli on a les relations suivantes : 

( £ } k = {^k}+ z ( K }k ( 7 - 16 ) 


Multipliant les deux membres de 1’ equation (7.16) par 




k} = 



({s°} + z{k}) 


(7.17) 


oil 




L’etat de contrainte dans le k e pli s’ecrit alors : 


1 

a 

I 

{e°} + z 

I 

a 

L J 

k 

_ _ 


(7.18) 


Les contraintes resultantes du stratifie peuvent etre obtenues a l’aide des relations entre 
les contraintes dans un pli et les contraintes resultantes : 


N 


<T 

X 

t 

X 

N v 

> = f v 


y 

J_i 

y 

X 

'Z 

2 

X 

L xy J 


t r_i t 

dz = j 2 r Q {s° jdz + j 2 t 


K 


}zdz 


En considerant que : 


(7.19) 
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[a 0 } et {k} du plan moyen geometrique sont independants de z: 

M 


JJ, 

Q 

|s°|dz = • 

v ) 

Jf, 

Q 

dz 

■{ e °} et \\ 

Q 

{ k} zdz = • 

A 

Q 

zdz 

2 


k 

2 


k 

V J * 

2 


k 

2 


k 


L Jk 


est constante suivant l’epaisseur de chacun des plis 


£ Q dz = Z Q f> q zdz= z Q f‘, zdz 

Jt k=l IL Jk J 2 L J k k=1 tL Jk 


Les integrates sont : 
f k dz 

Jz k-i 

1 


— Z k Z k-1 


I, 


Z( ^ Z ^ ( Z k Z k-l) 


(7.20) 


On obtient finalement 


N x 

< N y > 


Z- 

k=l 


( z k- z k-i) 


k=l IL Jk 


( z k - z k-i) w 




X' 


"A u 

A 12 

Ai6 


X 


X 

B 12 

B 16 


K x 

Ou 

« 

N y 

. = 

A u 

A 22 

A 26 

< 

s° y 

> + 

B 12 

B 22 

B 26 

< 

K y 



KJ 


_Ai6 

A 26 

A 66 _ 


y° 

L'xyJ 


_ B 16 

B 26 

B 66_ 


K *y. 


Avec 


N 

Ay = y.(Q„)k ( Z k 


k=l 


Z k-l) 


(7.21) 


(7.22) 


1 N 

B «=AZ(Q«) k ( z 

Z k=l 


k - z k-l) 


[A] est la matrice de rigidite en membrane et [B] est la matrice de couplage 
membrane-flexion-torsion. 
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De fagon similaire, pour les moments resultants, on a : 


[M x ] 


(7 

X 

t 

X 

M v 

> = r 2 ■ 

CJ„ 

y 

j_i 

y 

M *y. 

2 

t 

r-l 

X 

c 


t r_l t 

zdz = J 2 t Q js°]zdz + j 2 t 


2 L 


2 L 


K z dz 


M, 

m 5 

M 


xy 



L Jk 


(z 



<A 


J k-1 



soit : 


M x ' 


X 

B 12 

B 16 


X 


X 

D 12 

B 16 


K x 

M y 

> = 

B 12 

B 22 

B 26 

< 

B° y 

> + 

D 12 

d 22 

D 26 

< 

K y 

M xy 


_ B 16 

B 26 

B 66 _ 


r 

x o 

v 


_ D 16 

D 26 

D 66_ 


K xy. 


Avec : = 



1 N _ 

-ICQijMZk 

^ k=l 

1 N 

r£(Q,k(![ 

^ k=l 



(7.23) 


(7.24) 


[D] est la matrice de rigidite en flexion et torsion. 


Finalement, les relations entre les elements de reduction et les deformations du plan 
moyen geometrique peuvent s’ecrire : 


N x ' 


A n 

A 12 

■^■16 

B 11 

B 12 

B 16 


S ° 

N y 


A 12 

A 22 

-^26 

B 12 

B 22 

B 26 


By 

N xy 

. = 

•^■16 

A 26 

■^■66 

B 16 

B 26 

B 66 

< 

Yxy 

M x 

B 11 

B 12 

B 16 

Du 

D 12 

D 16 


K x 

M y 


B 12 

B 22 

B 26 

D 12 

d 22 

D 26 


K y 

X 

S 


_ B 16 

B 26 

B 66 

D 16 

D 26 

D 66 _ 


K xy. 


ou : 
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E 


' N ' 


" A 

B “ 


s 0 ' 


s 0 ' 

M 

> = 

B 

D 

« 

K 

> = E- 

K 










est la matrice de rigidite du stratifie 


(7.26) 


7.3 Matrice de rigidite des stratifies 

Stratifie symetrique 

Le pli a la hauteur z a les memes caracteristiques relatives a l'epaisseur, au materiau et a 
l'orientation que celui qui se trouve a -z. 

Remarque : La matrice [B] d’un stratifie symetrique est nulle. 


Exemple 7.1 Soit un stratifie [45/ -45/ -45/ 45] de composite carbone/epoxy AS/3501, 
dont les plis ont une epaisseur de t = 0.25 mm. Calculez la matrice de rigidite du 
stratifie. 


Solution 




= -0.50 
= -0.25 
= 0 

= 0.25 
= 0. 50 


Selon table (2.2) 


fEi =138 GPa 
[e 2 =9 GPa 


JG 12 = 6.9 GPa 
[v 12 = 0.3 


V 21 = v 12 = 0.3 x 


138 


0.0196 
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[Q] 


r e. 

v 12 E 2 

0 


"138.81 

2.72 

0 

1- v 12 v 21 

1- v 12 v 21 



E 2 

0 


2.72 

9.05 

0 


1 _ v 12 v 21 


SYM 


G 12 


0 

0 

6.9 


- V 12 v 21 = 1 - 0.3 x 0.0196 = 09941 


Qu = 


138 

1-0.3x0.0196 


138.81, 


Ql2 


q 22 - 

0.3x9 

0.9941 


9 

0.94941 
= 2.72 


9.48 


[q]= [T (6>)] _1 [q][t : (6>) 



“45.22 

31.42 

32.44“ 

GPa [Q]_ 45 = 

" 45.22 

31.42 

-32.44“ 

31.42 

45.22 

32.44 

31.42 

45.22 

-32.44 

32.44 

32.44 

35.6 


-32.44 

-32.44 

35.6 


A 11 - (Qll)l( z l _z o) + (Qh)2( Z 2 _z i) + (Qii) 3( Z 3 ~ z 2 ) + (Qll),4( z 4 _z 3) 


0.25 


0.25 


0.25 


0.25 


= 45.22 X (- 0.25 + 0.5) + 45.22(0 + 0.25) + 45.22(0.25 - 0) + 45.22(0.5 - 0.25) 
= 45.22[0.25 + 0.25 + 0.25 + 0.25] = 45.22 


[A] = 


“45.22 

31.42 

0 



“0 

0 

0“ 

31.42 

45.22 

0 

GPa - mm 

[B] = 

0 

0 

0 

0 

0 

35.6 



0 

0 

0 


GPa - mm" 


[d] = 


3.77 2.62 2.03 
2.62 3.77 2.03 
2.03 2.03 2.97 


GPa - mm ^ 
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Exemple 7.2 Prouvez que pour les stratifies symetriques J = 0 . 

Solution 


B ij=fi(Q.i) k ( z k- z k-i) 

z k=i 



Considerons le pli n° 1 et n 

Stratifie symetriquel z 0 = -z n 
pli n° 1 et n \ z i = -z n -i 


Symetrie => (q^ = (q^ = c te = C 

B U = K(°‘i) 1 ( Zl 2 " z “) + - + ( 0 l i)n( z ” 

d'oii : 

Bij =-|[(+ z n-i- z n)+- + ( z n- z n-i) 



= 0 


Stratifie equilibre / non symetrique (antisymetrique) 

Un stratifie equilibre et non symetrique est constitue de plis de meme epaisseur et qui 
sont disposes de fag on antisymetrique, ce qui veut dire que le pli au niveau z a une 
orientation +9 tandis que celle du pli au niveau -z est -9 . Notons que 6 ^ 0 et 0 ^ 90 . 

Exemple : [- 45 /- 45 / 45 / 45] 

Pour ce type de stratifie, on a les matrices de rigidite suivantes : 



'A u 

A 12 

0 “ 


" 0 

0 

B 16 


'Du 

D 12 

0 " 

[A]- 

A i2 

A 22 

0 

[b]= 

0 

0 

B 26 

[d] = 

D 12 

d 22 

0 


0 

0 

■^66 _ 


_ B 16 

B 26 

0 


^ 0 

0 

D 66_ 
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Exemple 7.3 Soit un stratifie [45 / — 45] de composite carbone/epoxy AS/3501dont les 
plis ont une epaisseur de 0.25 mm. Determiner la matrice de rigidite du stratifie. 


Solution 



-0.25 mm 
0 

0.25 mm 



"45.22 31.42 32.44“ 


45.22 31.42 -32.44“ 

[Q ] 45 - 

31.42 45.22 32.44 

GPa [Q ] -45 = 

31.42 45.22 -32.44 


32.44 32.44 35.6 


-32.44 -32.44 35.6 

A ij = Z(Qij ) k ( z k — z k-i)=(Qij ) 4J (t ) +(Qtj 

L(«) 


GPa 


k = 1 


ou t = 0.25 mm 

All = (o5 >45 (0 + (q5 )-45 ( t ) = 45.22 x 0.25 + 45.22 x 0,25 
= 22.61 GPa - mm = A 22 (a cause de Q n = Q 22 ) 
a 12 = ( 1 Q 12 ) 45 ( t ) + (Qi2 )-45 ( t ) = 31.42 X 0.25 + 31.42 x 0.25 


= 15.71 GPa - mm = A 21 

Am = A 26 = (O^ ) 45 (t ) +(o^)_ 45 (t ) = 32.44 x 0.25 - 32.44 x 0.25 

A 66 = (Q66 ) 45 (t ) + (Q66 )_ 45 ( t ) = 0 


= 35.6 x 0.25 + 35.6 x 0.25 = 17.8 GPa mm 



22.61 15.71 0 
15.71 22.61 0 
0 0 17.8 


N N 

GPa - mm - ou — x mm = 

mm 2 mm 
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Les coefficients de la matrice de couplage sont : = - £ ( Qj i ) k (z J - z\_ x 

2 k = i 

d’ou : 

i 


B 11 = 2 l Q ll )45 ( Z 1 - Z o) + (Qh )-45 ( Z 2 ~ Z 1 2 ) 

= - [45.22(0 - 0.25 2 ) + 45.22 (o.25 2 - 0)] = 0 = B 


B„ =- 


1 


[(Qi2) 45 fe-Zo)+(QnLfe-Zi 2 ). 


[31.42(0 - 0.25 2 ) + 31.42(o. 25 2 -o)] = 0 = B 2 


B 16 = j [32.44(0 - 0.25 2 )- 32.44 (o.25 2 - o)]= -2.027 GPamm 2 
= B, 


B 66 = ^ [( Q 6 6 ) 4 5 ( Z 1 " Z o) + (Q66 )_45 ( Z 2 - Z l 2 ). 


= ^ [35.6 x (o-0.25 2 ) + 35.6 (o.25 2 -0)] = 0 


[B] = 


0 0 -2.027 

0 0 -2.027 

2.027 -2.027 0 


GPa-mm 2 ou 10 3 N 


Rigidite en flexion : 


°U-| 

k=l 

d’oii 


ij )k V Z k Z k 




l z 2 - T>\ 


D 11 = ^ [( Q ll ) 4 5 ( Z 1 - Z o)+ (Qtt )-45 ( z 

45.22(0 - (-0. 25) 3 ) + 45.22 (o.25 3 - o)]= 0.471 


D 12 = 0.327 GPa - mm = D 2 i 

B 16 = B 26 = 0 
1 
3 


B 66 - 


T, .(.!() - 0.25 3 )+ 35.6(o.25 3 - 0 


= 0.371 GPa - mm' 
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et 


[D] 


0.471 0.327 0 
0.327 0.471 0 
0 0 0.371 


GPa - mm 3 


Stratifie quasi - isotrope 

Exemple [60/0/-60] et [90/45/0/- 45] 
Pour ce type de stratifie, on a : 


[A] 


1 

> 

A12 

0 

A i2 

A 22 

0 

0 

0 

An — A 12 

2 




Exemple 7.4 Prouvez que le stratifie (0/90) n’est pas necessairement quasi-isotrope. 

Solution 




Qn 

Q12 


o 


Ql2 0 

q 22 0 



0 



Q 22 q 12 0 

Ql2 Qn 0 
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A ij =S(Qij) k ( z k - z k-i) 

k=l v ' 

t 

A 11 = (Qh ) 0 (0 + (Qh ) 90 (0 = Qll t + Q22 t = (Qn + Q22) t 

Supposons que le stratifie a une epaisseur t = 1 , on obtient : 

A u = (Q n + Q 22 ) 

A 12 = (Ql2 ) 0 (0 + (Ql2 ) 90 (0 = 2 Ql2 t= 2 Ql2 
A 6 =0 

A 22 = (Q 22 ) 0 (0 + ( Q 22 ) 90 (0 = (Q 22 + Qll) t = (Q 22 + Qll) 


0 

0 

:Q66_ 

(1) 

( 2 ) 

(3) Ok 


(4) pas necessairement 

1 Ml 1 ^12 

66 _ r— » 


A 66 - (Qb6 ) 1 + (Qb 6 ) 90 1 “ 2 Qe6 

(Qi 1 +Q 22 ) 2 Ql 2 

[ A ]= 2Q 12 (Q11+Q22) 

0 0 

Pour un stratifie quasi-isotrope, on a 


A 11 A 22 


An — A 17 
et A cc = ^ 


La matrice [a] donne : 
A 11 = (Qn + Q22) = A 2 


a = ?n Q11 + Q22 ~ 2 Qn 
-^66 Z ^66 
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7.4 Matrice de souplesse des stratifies 

Inversion de la matrice de rigidite 

Les relations entre les forces resultantes et les deformations sont : 



' N ' 


" A 

B “ 


[e°] 


fe°] 


< 

M 

" — 

B 

D 

< 

K 

= [E] 

K 

> 


Inversement 



(7.27) 


(7.28) 


La matrice de souplesse d’un stratifie est l’inverse de sa matrice de rigidite. On peut 
inverser la matrice de rigidite [E] du stratifie soit directement soit par une methode 
alternative (par partition) demontree par Halpin : 


[E] 


“A' 

B' 

B' 

D' 


oil : 


[a'Ha-Hb-Id11c-] 

[B'Hb'IdT 

[c] = - [d’ ] 1 [c' ] = [b 1 ] 
[d']=[d’1‘ 


et 

A* 

B* 

c* 

D* 


[A] 1 

-[ArM 

[BlA ] 1 

[D]-[BlAnB] 


(7.29) 
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Calcul des constantes elastiques equivalentes de stratifies 

Les constantes elastiques sont parfois plus utiles que les matrices de rigidite [A], [B] et 
[D]. Ces constantes peuvent se determiner en utilisant les matrices [A 5 ], [B’] et [D 5 ]. 

Stratifie symetrique soumis au chargement dans son plan 

Dans ce cas, [B] = 0, d’oii{^ x } = {^ 0 }. Par consequent, le champ des deformations 
devient : 


s x 


r 

O X 
CO 

v. 


£ y 

> = < 

B° y 

. = 

.V 


Y° 

I'xyJ 



An -A-12 -A- 16 

A ' A ’ A ’ 

-^12 -^22 -^26 

A ' A ’ A ’ 

-^16 -^26 -^66 



f N ] 


X 

< 

N 


y 


>> 

X 

£ 

J 


(7.30) 




Figure 7.7 Stratifies soumis aux chargements dans son plan [Gibson] 


Chargement 


N x 

0 

0 


=> = A u N x , * v ° = A i2 N x et Xx ° y = A i6 N 


16 x 


Les constantes elastiques sont : 
N„ 


F 


xy 


,° 

x y 


t 

1 

(7.31) 

A n N x 

tA n 

_ A i2 

An 


(7.32) 

_ A i6 
"A n 


(7.33) 


Chargement 


0 

N„ 


' => £ y = A 22 N y 
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N, 


E = G ^- = 


y e° y A 22 tA 


(7.34) 


22 


Chargement 


0 

0 


N 


xy 


> £°v = A 26 N xy et y° v = A 26 N 


xy 26 xy 


N 


xy 


'xy 


° v = y% A 


t 


tA c 


(7.35) 


Tlxv.v 


£y _ A 2 g 


'xy,y o 
xy 


y xy -^-66 


(7.36) 


Stratifie symetrique soumis a un moment de flexion 

Les relations entre les rayons de courbure et les moments sont : 


K x 


1 

o_ 

d; 2 

D 16 


X' 

s 

. = 

D12 

D 22 

^26 

< 

M y 

K xy 


•^16 

^26 

^66 


M xy. 


(7.37) 



Figure 7.8 Stratifies en flexion pure [Gibson] 
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car k = 


Chargement < 
1 


M x 

0 

0 


. „ , M 1 

1=>k x = D u M x =D U — = — 

b Px 


or 


Px 

1 M 


M 


Px Efx^yy g bt^ 

fx 12 

Le module effectif en flexion du stratifie suivant la direction x est done 

E f = 

t 3 D„ 


(7.38) 


De la meme faqon, on obtient : 
12 


E t . =■ 


t D. 


(7.39) 


Verification experimentale 

La theorie classique des stratifies peut etre verifiee experimentalement en mesurant les 
deformations dans un stratifie qui est soumis aux chargements connus. Ces valeurs 
mesurees seront par la suite comparees avec les valeurs calculees a l’aide de la theorie. 
Tsai a propose cette methodologie afin de verifier les resultats theoriques des stratifies de 
composite verre-epoxy [Tsai, 1964]. 


Determination des deformations et courbures 
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Collons une rosette a la surface superieure et une autre sur la surface inferieure de 
l’eprouvette. 

A la surface superieure, les deformations suivant les directions x et y sont : 



(7.40) 


(7.41) 

A la surface inferieure, on a : 


sf =e“ +^k x 

(7.42) 

e“ = 8° + — K 

y y 2 y 

(7.43) 

(7.40) et (7.42) => s° et k x 

(7.41) et (7.43) => ej et K y 


La deformation suivant la direction 0 s’ecrit comme suit : 


s s e up = s s x up cos 2 (0) + s s y up sin 2 (0) + y x up cos (0) sin(0) 

(7.44) 

pour 0 =45° la deformation angulaire est : 


(y7) 45 = 2s7-(8 s x up +s s ; p ) 

(7.45) 

( Y xy ) 4 5 = 2s^5 - (£x f + e y f ) 

(7.46) 

or : 


y sup =y° — — k 

1 xy lx y 2 xy 

(7.47) 

y inf = y° +1 K 

1 xy l xy 2 xy 

(7.48) 

La resolution du systeme comportant ces dernieres 4 equations, on obtiendra : y° y et K xy . 

En utilisant cette eprouvette instrumentee des jauges de contraintes, il est done, possible 
de determiner les coefficients des matrices [A’], [B’] et [D’] a partir des donnees 
experimentales. 
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Determination des matrices [A’] et [B’ ] a Vaide d’un essai de traction 

Dans un essai de traction suivant la direction x, les deformations du PMG d’un stratifie 
sont : 


'4 


A n 

A i2 

1 

CD 

- 1 

< 


N x ' 

O 

CO 

> = 

A i2 

a 22 

A 2 6 

« 

0 

Y° 

L'xyJ 


Ai6 

A 2 6 

A66 


0 


4=AuN x => A n =-i 

^ X 
0 


s° y =A 12 N, 


=> A 12 - 


N, 


r" y =A; 6 N x => a; 6 =^l 
Les courbures causees par N x sont : 


K x 


Bn 

B i2 

Bl6 


N x ' 

K y 

. = 

B i2 

b 22 

B 2 6 


0 

K xy. 


Bi6 

B 2 6 

B66 


0 


(7.49) 


K x = B n N x 
K y = B i2 N x 
K y = Bi2 N x 


=> B ii = 




N x 


- > B 12 - 


K, 




_ > B 16 - 


K 


xy 




(7.50) 


Determination de [D’] a Vaide d’un essai de flexion pure 


Chargement : 


M x 

0 

0 
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Les coefficients de la matrice [D’] sont : 

■ k y 

D n= — 

M y 


D12 - 


K, 




K xv 

n — y 
u 16 ~ 




(7.51) 


7.5 Calcul des contraintes et deformations dans les plis 

Dans le k e pli, les contraintes sont : 


{ a (*,y)} k - 

Q 

({ £ (x,y)} + z k{ K (x,y)}) 



k 



[e°] 


" A 

B “ 

-l 

' N ' 


“A 1 

1 B' 


' N' 

ou • 

K 

" 

B 

D 

< 

M 

* — 

B' 

1 D' 


M 


(7.52) 


Exemple 7.5 Un stratifie [45/ -45] de composite carbone/epoxy AS/3501 dont les plis 
ont une epaisseur de 0.25 mm est soumis a une force de traction uniaxiale N x = 30 MPa 
- mm ou N/mm, calculez les deformations et les contraintes dans les plis du stratifie. 

Solution 

A partir des resultats de 1’ exemple 7.8, on obtient la matrice de rigidite : 




"22.61 

15.71 

0 

0 

0 

-2.027 " 



15.71 

22.61 

0 

0 

0 

-2.027 

“A 

UB" 


0 

0 

17.8 

-2.027 

-2.027 

0 

B 

1 D 


0 

0 

-2.027 

0.471 

0.327 

0 



0 

0 

-2.027 

0.327 

0.471 

0 



-2.027 

-2.027 

0 

0 

0 

0.371 


Unites = Ay (GPa - mm) 
Bij (GPa -mm 2 ) 
Dy (GPa -mm 3 ) 
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1 


GPa - mm 


: (MPa - mm)x 10 3 GPa 


MPa 


' 0 



£ x 


30 

/ 

0.003043 mm /mm 

£ y 


0 


-0.001183 mm /mm 

x o 

▼ 

•=[Er‘- 

0 

0 

■xl0“ 3 =• 

0 

0 

K x 





0 


0 

K y 


0 


0.01017mm” 1 

[ K xyJ 




Pli 

Location (mm) 

a, (MPa) 

a y (MPa) 

G y (MPa) 

N° 1 sup. 

-0.25 

18.74 

-41.26 

-30.18 

N° 1 inf. 

0.0 

101.3 

41.26 

60.36 

N° 2 sup. 

0.0 

101.3 

41.26 

-60.36 

N° 2 inf. 

0.25 

18.74 

-41.26 

30.18 




-PI 


K=1 

Superieur 


K=1 

Superieur 


avec 


£ x 




ZK x 


£ y 

> = < 

£ y 

> + • 

ZK y 

> 

y xy 

Superieur 

Yxy. 


ZK xy. 

Superieur 
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7.6 Influence de l’humidite et de la temperature 

Degradation des proprietes mecaniques 

L’humidite et la temperature affaiblissent les proprietes des composites notamment Ei, 
E 2 , G 12 et vi 2 . Par consequent, les proprietes mecaniques d’un stratifie qui est expose a 
un environnement chaud et humide seront affectees par la chaleur et par l’humidite. 


Contraintes hygrothermiques 

Dans le k e pli, les deformations totales comprennent les deformations mecaniques, 
hygroscopiques et thermiques : 



(7.53) 

par consequent, les contraintes sont : 



(7.54) 

Dans le cas d’un pli, lorsque les deformations totales {s} = 0, les contraintes 
hygrothermiques sont : 


CT = 


(-{<x}AT-{p}C) 


(7.55) 


Pour les stratifies, en general {s} k = {s° }+ z{k} ^ 0 

Les contraintes totales incluant les contraintes hygrothermiques deviennent alors : 



({ £ °} +z { K } - {a} k AT+(p} k C) 


(7.56) 


Les contraintes resultantes totales peuvent s’ecrire : 
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{ N (x,y)} = / 2 thx,y)} k dz 


f 2 1 Q ({ £ (x,y)} + z { K (x,y)} { a (x,y)} k AT {P(x,y)} k C ) dz 


2 J k 
0 


-\A\\ ej'x.y) | + [B J ji^fx.y) } - { N (x,y) } “ { N (i,y) } 


oil les contraintes resultantes d’origine thermique sont : 


N 


( x .y) 




t r _ n 

Q 


a 


(x.y) 


ATdz = AT^ 


2 L J k 


k=lL Jk 


a (x,y)ji,( z k ~ z k-l) 


et les contraintes resultantes hygroscopiques sont : 


N " ,y)} = \\ [«] k {P(x,y)} k Cdz = C g [q] {P(x,y) } k < Z K “ z K-i> 


k=lL Jk 


De fagon similaire, les moments resultants peuvent s’ecrire : 

' '•» v 


{ M (*,y)} = / 2 tK.y)} k zdz 


t 


j 2 1 ({ e (x,y)} + z { K (x,y)} {“(x,y)} k AT {P(x,y)} k c ) zdz 


[ B ]{ e (x,y) | + [D]K X ,y)}- 1 <y) } ~ Kx,y) 1 


oil les moments resultants thermiques sont : 


t r ~\ 


M 


( x »y) 


f 2 , 


a 


(x.y) 


2 L J k 


ATzdz = ^Z[Q] k { a (x,y)} k ( Z k ~ Z iU) 
z k-1 


et les moments resultants hygroscopiques sont : 


M 


H 

(x.y) 


2 L -Ik 


i-1 Q _ {P(x,y)} k Czdz “ 2 ij[Q]k{P(x,y)} k ( Z k Z k-l) 


H 2014 


(7.57) 


(7.58) 


(7.59) 


(7.60) 


(7.61) 


(7.62) 
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En rearrangeant les equations (7.57) et (7.60), on obtient les relations suivantes 


{ N 0.y>! + { N (x,y> } + - [ A ] { S °K.y)} + [ B ]{ K 0.y)} 

K,,} + K,») + W,,} = [ B ]{4.y,} + [°]{ Vy>) 


ou encore : 




" A 

B “ 

<| 

fo° 1 
£ (x,y) 

M fx,y) 


B 

D 


\ 

>< 

j 


oil les contraintes et les moments resultants effectifs sont 

K«)J = { N o.y)} + H.y)i + { N ".,)} 


M (ky)} “ { M «,y>} + ( M 0<y)} + { M ('x,yl} 


(7.63) 


(7.64) 


(7.65) 


Remarque : 

{N} et {M} sont les contraintes et les moments resultants mecaniques externes. 


Les deformations sont obtenues en inversant la relation (7.64) : 


fo0 1 

£ (x,y) 


“ A 

B “ 

-l 

N(x,y) 


"A' 

B' 


N (x,y)' 

K (x,y). 


B 

D 


M fx,y). 


B' 

D' 

< 

M fx,y), 


(7.66) 


Exemple 7.6 Un stratifie [-45/ 45/ -45/ 45] de composite carbone/epoxy AS/3501 dont les 
plis ont une epaisseur de 0.25 mm, est chauffe de T 1 = 20 °C a T 2 = 100°C. Supposons 
que les proprietes ne sont pas influences par la temperature, calculez les contraintes 
thermiques engendrees. 

Solution 



' Zq = — 0. 50 
z, = -0.25 


^ =0.25 
z 4 = 0. 50 
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Tab. 5.2 => 


oq = 0.88 xlO^ 6 /°C 
< a 2 =31.0x10 6 /°C 
a 12 = 0 


xy 


> = [T , (0)] 1 . 

oq 

a 2 


, a i 2 . 


a x 

a y 

, — 

T (45) 

-1 

< 

0.88 

31.0 

>xlO“ 6 /°C = - 

15.94 

15.94 

Cl 

^xy 

45 



0 


30.12 


>xlO“ 6 /°C 


a. 


a. 


a 


N 


xy 
T 

( x .y) 


— 

T (-45) 

-1 

« 

0.88 

31.0 

>xlO“ 6 /°C = - 

15.94 

15.94 

/IC 



0 


-30.12 


>xlO“ 6 /°C 


-45 


N 


x2x0.25 = < 


AT Z[Q] k { a (x,y)} k ( z k-z k -i) 

k=l 

(x,y)} = ( AT )([Q]„45 { a (x,y)}_ 45 + [Q] 45 { a (x,y) 

: ~Y ( [^]-45 ( a ( x -y) } 45 (Z 1 “ Z 0 ) + PL { a ( x >y) } 45 (Z 2 “ Z 1 ) 

+ P-45 ( a ( x >y) }— 45 (Z 3 “ Z 2 ) + [Q] 4S { a (x,y) } 45 ( Z 4 “ Z 3 ) = 


19.5755 

19.5755 

0 


> MPa - mm 


0 

0 

-0.3834 


MPa — mnT 


0.000307 

0.000307 



0 

-0.000967 
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Contraintes : 


:a }k = [Q] k ({ s °} +z W _ Wk AT 



-11.764 



-3.921 ' 

{a} =< 

i ) lsup eneure 

-11.764 

-5.739 

MPa 

^Jlinf erieure 

-3.921 

-14.348 


>MPa 


,etc. 


7.7 Resistance du stratifie 

Resistance du pli le plus faible 

Remarque : Les contraintes inter-plis ne sont pas considerees. 



Figure 7.10 Courbe « contrainte-deformation » d’un stratifie soumis a une charge 
de traction indiquant des ruptures successives des plis conduisant a sa rupture 
finale [Gibson] 
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Pour la n e section de la courbe « contrainte-deformation », on a les relations suivantes : 


N &) 


A (n) 

B (n) 

<i 

f P ° (n) 1 
8 (x,y) 



B (n) 

D (n) 


k- 

[ K (x,y) J 


(7.67) 

Les elements de reduction totale sont : 


N (x,y) 

> 

ii 

N 

iN (x,y) 

M (x,y) 

tot 

n=l 

[ M (x,y) J 


(7.68) 

et les deformations totales : 


.. 


fp° (n)l 

8 

k 

£ (x,y) 

— 

II 

M 

— 

K 

n=l 

K (n) 

l K (x,y) J 


tot 


7.8 Theorie des plaques stratifiees du premier degre (First Order Shear 
Deformation Theory) 


Deformations 

De faqon generale, les champs de contraintes et de deformations dans une plaque mince 
constitute de stratifies de composites peuvent se decrire avec une bonne precision en 
utilisant la theorie classique des stratifies (TCS), a l’exception des regions pres de leurs 
bords libres. Par contre, pour la description du comportement mecanique (fleche, 
repartition des contraintes, etc..) des stratifies epais dont le rapport entre la largeur et 
l’epaisseur est inferieur a 10, cette theorie devient mal adaptee. La theorie des plaques du 
premier degre est une amelioration en tenant compte du cisaillement transverse . 

Hypotheses : 

La deformee d’une normale AA’ du stratifie reste un segment de droite lors de la 
deformation. Toutefois, contrairement a la theorie classique des stratifies, cette deforme ne 
reste pas normale a la deformee du plan moyen geometrique (PMG) (figure 9.2). 
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La longueur de AA’ reste inchangee; ce qui implique que s z =0. 



Figure 7.11 Deformation selon la theorie des plaques du premier degre en tenant compte 

du cisaillement transverse. 


Sous l’action de la charge externe, le point p° qui est situe dans le PMG se deplace d’une 
distance u° suivant la direction x (figure 7.11). Le point p qui est situe d’une profondeur z 
en dessous du point p° se deplace d’une distance : 


u(x, y, z) = u° (x, y) - z(j) x (x, y) 

(7.70) 

v(x,y,z) = v (x, y) -zcj) y (x, y) 

En considerant l’hypothese (2) et le fait que les angles cp( X ) et cp( y ) sont tres petits, il est 
admis que : 


w(x, y,z) = w°(x,y) (7.71) 

En utilisant les relations entre deformations et deplacements, les composantes de la 
deformation en un point p quelconque de l’epaisseur ont pour expression : 
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du 8 o d(j) du° d(j) x o 

£ * = ^r = ^ (u ~ z ^f )= ^r~ z ^ =£ * +ZK * 

dx dx ox ox ox 


dv 8 o S(j) dv° d(j) 

5 = = (v -z — -) = z — 

dy dy dx dy dy 


£ y +ZK y 


du dv du° dv° dcj) x dcj) y 


Yxy = ^T + ^ = ^r + ^ Z (^ 2L + ^ JL ) = Yxy+ZK xy 

dy dx dy dx dy dx 

dv dw , dw° 

y — 1 — — (b H 

^ dz dy y dy 

du dw , dw° 

Yxz + ~~Vx — 

dz dx dx 


(7.72) 


ou : 


5( K 

dx 
dy 

K = _A + ^ il 

xy dy dx 


K = 


K y=“ 


(7.73) 


Lorsque l’epaisseur de la plaque est tres faible par rapport aux autres dimensions, les 
deformations y yz et y xz peuvent etre negligees : 


. dw° n . dw° 
1yz y dy y dy 

_ . dw° _ , _dw° 

dx dx 


(7.74) 


En substituant (7.74) dans (7.73), les expressions des courbures deviennent identique aux 
celles dans le cas de la theorie classique des stratifies : 


a 2 o 
d w 

Kx_ dx 2 

a 2 o 
d w 

Ky_ ~d?~ 

d 2 w° 
xy dxdy 


(7.75) 
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Contraintes dans le k e pli 

Contraintes en membrane : 


<7 




8 V 

X 




X 

a y 

. = 

Q 

< 

£ y 




k 


.V 

k 



Yxy, 


Contraintes du cisaillement transverse : 



(7.76) 


(7.77) 


Matrices de rigidite et de souplesse du stratifie 



Figure 7.12 Les contraintes et moment resultantes 

Les contraintes resultantes en membrane et moments resultants sont definis dans la 
section 7.2 comme suivantes: 
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N xy =j 2 t T xy dz-^|j 2 t (T xy ) k dz| 

2 k=l [ 2 J 

M x = f 2 t^x Zdz = Z|f 2 tK)k Xzdz } 

2 k=l l 2 J 

M y = j " °y ZdZ = X ( fe K )k X Z dZ | 

(T xy )k Xzdz 



(7.78) 


d’oii : 


X' 


A n 

A i2 

■^■16 

B 11 

B 12 

tn 

1 


U] 

N y 


A i2 

A 22 

■^■26 

B 12 

B 22 

B 26 


< 

N xy 


A-16 

■^■26 

■^■66 

B 16 

B 26 

B 66 


Yxy 

> — 







< 

M x 


B u 

B 12 

B 16 

D n 

D 12 

D 16 


K x 

My 


B 12 

B 22 

B 26 

D 12 

d 22 

D 26 


K y 

KJ 


_ B 16 

B 26 

B 66 

D 16 

D 26 

D 66_ 


K xy. 


avec : 

N 

Ajj = y.(Qii) k ( z k z k~i ) 

k=l 

B iJ=IZ(<W z G z L) ( ? 80 ) 

^ k=l 

■J k=l 


Par analogie, les contraintes resultantes du cisaillement transverse presentees a la figure 
7.12 sont definies comme : 
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(7.81) 


d’oii : 



(7.82) 


avec : 


F a =Z(Cfl) k ( z k-Vi) iJ = 4 - 5 


(7.83) 


Remarque : En general, les matrices [A]. [B]et [D] sont utilisees dans l’application de la 
theorie classique des stratifies. La matrice [F] qui s’appelle matrice de rigidite en 
cisaillement transverse, est utilisee seulement lors de l’application de la theorie des 
plaques de premier degre. 

7.9 Theorie modifiee des plaques stratifiees avec cisaillement 
transverse 

II s’agit d’un remplacement des coefficients de la matrice de rigidite des contraintes du 
cisaillement transverse [F] par de nouveaux coefficients. La nouvelle matrice resultante 
est la matrice [H] : 



(7.84) 


avec 



(7.85) 


Les parametres ky sont les facteurs de correction en cisaillement a determiner. En 
supposant que la repartition de la contrainte de cisaillement est une parabole en fonction 
de l’epaisseur du stratifie, la fonction suivante peut etre utilisee pour la ponderation : 


f(z) = (5/4)[l-(2z/t) 2 ] 


(7.86) 
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Par consequent, la matrice H s’ecrit : 

H .)=TE( e «) k [(l'i-h k _ 1 )-i(h=-h 3 l _,)i] (7.87) 

4 k=l 3 t 

oil t = epaisseur du stratifie et i,j = 4,5. 

Exemple 7.8 Soit un stratifie [55/-55J forme de deux plis de materiau composite 
unidirectiormel dont C44=G 2 3=0.385GPa et C 5 5=G 13 =3.789GPa. Sachant que les plis 
out la meme epaisseur de 0.635mm, calculez la matrice de rigidite reduite pour les 
contraintes du cisaillement transverse. 


Solution 

Pli #1 (9=55 -1 

(C 44 ) 55 = C 44 cos 2 0 + C 55 sin 2 0 = 0.385cos 2 (55) + 3.789sin 2 (55) = 2.669GPa 
(C 55 ) 55 = 0.385 sin 2 (55) + 3.789 cos 2 (55) =1.505GPa 
(C 45 ) 55 = (3.789 - 0.385) sin(55) cos(55) = 1.599GPa 

Pli #2 (9=-55 -) 


(C 44 )_ 55 = C 44 cos 2 0 + C 55 sin 2 0 = 0.385 cos 2 (-55) + 3.789 sin 2 (-55) = 2.669GPa 
(C 55 )_ 55 =0.385sin 2 (-55) + 3.789cos 2 (-55)=1.505GPa 
(C 45 )_ 55 = (3.789 - 0.385) sin(-55) cos(-55) = -1.599GPa 


En appliquant l’equation (7.87), on obtient : 

H 44 = 2.82GPa-mm 
H 55 = 1.59 GPa-mm 

h 45 = o 
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Exercices 

1. Un stratifie dont P epaisseur est egale a 0,25 mm, est soumis a un chargement dans 
son plan. Sachant que les deformations du plan moyen geometrique sont les 
suivantes : 


A 


2751 




0,077 


< 

■ — - 

-1331 

xlO' 6 et 

K °y 

- = ■ 

0,009 

mm 1 

lly 


-1125 




0,070 



Determinez les deformations a la surface superieure du stratifie. 


Reponse : 


0 

£ x 


6874 

0 

£ 

- = - 

-2456 

y 



0 

Ay. 


-7625 


xlO 


2. La figure suivante illustre un stratifie [0/30/-45] de carbone/epoxy AS3501 (vf = 
65%). Sachant que les plis ont une meme epaisseur de 5mm, calculez la matrice de 
rigidite du stratifie. 



1,34 

0, 29 

0,04 

Reponse : [a] = 

0, 29 

0, 37 

-0, 08 


0,04 

-0,08 

0,35 


xl0 6 MPa-mm 



-2, 34 

0,72 

-0,81 



2,58 

4,87 

-3,97 

B] = 

0, 72 

0, 90 

-0,81 

X 10 6 MPa-mm 2 , 


4,87 

7, 56 

-4,23 


-0,81 

-0,81 

0, 72 



-3,97 

-4,23 

6,05 


3. Soit un stratifie [±30/0] s de carbone/epoxy dont les caracteristiques sont les 
suivantes : 


Ei = 150 GPa, E 2 = 12,1 GPa, G 12 = 4,4 GPa, v 12 = 0,248. 

Sachant que les deformations du plan moyen geometrique sont : 

6° = 1000 xl0~ 6 ete y = = k,° = k° — K° y = 0 , determinez : 

les chargements causant ces deformations; 

la repartition des contraintes a travers P epaisseur du stratifie. 

Reponses : N x = 99,2164 N/mm, N y = 18,828 N/mm, N xy = M x = M y = M xy = 0; 

Pli#l,sup. : ct x = 89,9866(MPa), cr y = 29,1301(MPa), t xv = 29,1301(MPa) 

Soit un stratifie [0/90] s d’un composite carbone/epoxy dont les 
caracteristiques fondamentales sont les suivantes : 
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Ei = 181 GPa, E 2 = 10,3 GPa, G 12 = 7,17 GPa, v 12 = 0,28 
Sachant que l’epaisseur de chacun des plis est de 5 mm, calculez : 

E x , E y , G xy , v xy , Ef X et E fy du stratifie. 

Reponses : E x =124,53 GPa, E y =67,43 GPa, G xy =7,17 GPa, v xy =0,043, Ef X = 174,96 GPa, Ef y =16,65 GPA 
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